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Precise cks from r Decays 

B. Malaescu'" * 

'^Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11 (UMR 8607), 
F-91405, Orsay Cedex, France 

An updated measurement of as(m^) from ALEPH r hadronic spectral functions is presented. We report a 
study of the perturbative prediction(s) showing that the fixed-order perturbation theory manifests convergence 
or principle problems not presented in the contour-improved calculation. Potential systematic effects from quark- 
hadron duality violations are estimated to be within the quoted systematic errors. The fit result is as{m^) — 
0.344± 0.005 ±0.007, where the first error is experimental and the second theoretical. After evolution, the as{m%) 
determined from r data is one of the most precise to date, in agreement with the corresponding N^LO value 
derived from Z decays. 



1. INTRODUCTION 

The hadronic decays of the r lepton provide 
a clean laboratory to perform precise studies of 
QCD. Invariant mass distributions obtained from 
long distance hadron data allow one to compute 
the spectral functions, which permit the study of 
short distance quark interactions. In particular, 
these spectral functions can be exploited to pre- 
cisely determine the strong coupling constant at 
the r-mass scale, as(m^) . Most of the present 
analysis is described in detail in ref [1] . Some new 
studies on the perturbative methods are described 
in Sec. [H 



2. TAU HADRONIC DATA AND SPEC- 
TRAL FUNCTIONS 

The nonstrange vector (axial-vector) spectral 
functions ui(ai), for a spin 1 hadronic system, 
are obtained from the squared hadronic mass 
distribution, normalised to the hadronic branch- 
ing fraction (with the hadronic observable (s) 
6 - - 

Rt- v/a — ^'^ '^'"^ ), and divided by a factor 
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exhibiting kinematics and spin characteristics 
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The basis for comparing a theoretical description 
of strong interaction with hadronic data is pro- 
vided by the optical theorem, which relates the 
imaginary part of the polarisation functions on 
the branch cut along the real axis, to the spectral 
functions: 2-k ■ /TOn^^^(s) = Vijaxis). 

The determinations of i?T- from measured lep- 
tonic branching ratios, or only from the elec- 
tronic one assuming universality, are in very good 
agreement, yielding Rt = {1 — Be — Bf_i)/Bc = 
l/i3^"'- 1.9726 = 3.640 ±0.010. One can identify 
in Rr a component with net strangeness and two 
nonstrange vector(V) and axial-vector(A) compo- 
nents. Including the latest results from BABAR 
and Belle the value of the strange component is 
R^ S = 0.1615 ± 0.0040 . The separation of the V 
and A components for final states with only pi- 
ous is done using G-parity. For the KKtt mode, 
we assume CVC and use new results from the 
BABAR Collaboration [2j. Finally, we get the 
components: Rr,v = 1-783 ± 0.011 ± 0.002 and 
Rr,A = 1-695 ± 0.011 ± 0.002 , where the first er- 
rors are experimental and the second due to the 
V/A separation. 
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3. THEORETICAL PREDICTION OF 

The nonstrange ratio Rry/A can be written as 
an integral of the spectral functions over the in- 
variant mass-squared s of the final state hadrons 



RT,v/A(sa) oc / — ^mWyj^is + ie) 



(2) 
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The two point correlator can not be predicted by 
QCD in this region of the real axis. However, 
using Cauchy's theorem, one can relate this ex- 
pression to an integral on a circle in the complex 
plane. Then, the OPE yields Rr,v/A « 1 + d^"'> + 

^EW + ^udy'/A + Ed=4,6,... ^i%/A ' witli a mass- 
less perturbative contribution, a non-logarithmic 
electroweak correction, the dimension two per- 
turbative quark-mass contribution and higher di- 
mension nonperturbative condensates contribu- 
tions respectively. The perturbative part reads 
5(0) = Knk)A^"Has) , with the fimctions 



cause important problems as the absolute value 
of T] gets large and the convergence speed of the 
series is reduced [T]. In addition, a cut at a fixed 
order in as(so) is applied on the Taylor series and 
on the integration result in FOPT. Therefore, im- 
portant known higher order terms are neglected, 
yielding additional systematic uncertainties. 

A better suited method is CIPT which, at each 
integration step, computes as{s) using the value 
found at the previous step. In this approach the 
Taylor expansion is always used for small abso- 
lute values of its parameter, hence excellent con- 
vergence properties. 

One can analytically prove that the FOPT re- 
sult can also be obtained by making small steps, 
with a fixed order cut of the result at each step. 
However, in that case the effective RGE is mod- 
ified at every single step. This is another way 
to see the problems of the FOPT method, which 
exist for a transformation on the circle in the com- 
plex plane, as well as for a scale transformation 
on the real axis. 

At the order /3o one can analytically compute 
the solution of the RGE, 
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where Ug = as/ir and ^ is a scale factor. A 
breakthrough was made recently [5], so that the 
pertubative coefficients are now known up to 
K4 (see pP for the numerical values of the Kn{(,) 
coefficients) . 

3.1. Perturbative methods 

The perturbative contribution to Rr provides 
the main source of sensitivity to as (so). The 
value of the strong coupling in the complex plane 
can be computed assuming the validity of the 
renormalisation group equation (RGE) outside 
the real axis, and using a Taylor series of 77 = 
ln{s/so). 

In the fixed order perturbation theory (FOPT), 
at each integration step, the Taylor expansion is 
made around the physical value as(so). This can 
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On the real axis, as well as on the circle in the 
complex plane, CIPT reproduces the analytical 
solution 11]), whereas FOPT its approximation ([5]). 

At higher orders, one can analytically compute 
the integral of the inverse beta function, yielding 
at the order (3i 
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and then relate it to ln{s/so) using the RGE. 
Solving the resulting equation one gets a solu- 
tion that we will call af^{s). It is also instructive 
to compare the results of the two methods with 
a solution proposed for example in ^ , consisting 
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of an expansion in inverse powers of logarithms 
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Figure 1. The evolution of as{s) to lower scales 
^TO^ using CIPT, FOPT, the analytical solution 
of the RGE {a'^"{s)) and the expansion in inverse 
powers of logarithms {a^^'~^ {s)). These solutions 
were computed up to the order f3i of the RGE 
(top), and /32 (bottom). 

In Fig. [T] we show a comparison of scale trans- 
formations performed with CIPT, FOPT and the 
two other RGE solutions described above. CIPT 
yields to the same result as the analytical solu- 



tion, and the expansion in inverse powers of log- 
arithms is closer to it than to FOPT. We should 
emphasize the fact that, on the contrary to CIPT, 
the FOPT solution is problematic even on the real 
axis, generally not satisfying to the RGE. 

For numerical applications we have used geo- 
metric growth estimations for the first unknown 
coefficients /34, and Kq. We have tested that 
CIPT is less sensitive to changes of these coef- 
ficients, and it also exhibits a smaller scale de- 
pendence than FOPT. Numerically, the difference 
of the perturbative contributions computed with 
the two methods is about 15%. In fact this dif- 
ference could have been much larger if not for the 
properties of the kernel in the integral ([3]) which 
has small absolute values in the region where the 
as{s) predictions of the two methods are rather 
different [1]. 

In ref. [1], we analysed further problems with 
the perturbative series for FOPT when the in- 
tegration over the circular contour is performed. 
CIPT behaves better than FOPT and is to be 
preferred. The difference between the results ob- 
tained with the two approaches is not to be in- 
terpreted as a systematic theoretical error, but 
rather like a problem of FOPT [1 i 

3.2. Quark-hadron duality violation 

It is known that OPE describes only part of 
the nonperturbative effects. In order to estimate 
the impact of the missing contributions, we test 
two models based on resonances and on instan- 
tons. We add their contributions to the theoreti- 
cal prediction, choosing parameters that provide 
a good fit to the V+A spectral function near the 
r mass. For these models, we find corrections 
situated within our systematic uncertainties [1]. 

4. COMBINED FIT 

In order to obtain additional experimental in- 
formation, we use spectral moments defined as 







(8) 



^ See however a recent study with different conclusions [5] . 
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They allow one to better exploit the shape of the 
spectral functions and they suppress the region 
where OPE fails. The corresponding theoreti- 
cal prediction is very similar to the one for Rr , 
with consequent perturbative and nonperturba- 
tive contributions. Due to strong correlations, we 
use only Rr {k = and I = 0) and four additional 
moments (fc = 1 and I = 0,1,2,3) to simulta- 
neously fit asirn^) and the leading D = 4, 6, 8 
nonperturbative contributions. 

In spite of the fact that the nonperturbative 
contributions fitted for the V and A spectral func- 
tions have opposite signs and they are much larger 
than those from V+A, we find an excellent agree- 
ment between the values found for Q;s(m^) from 
the three fits. The result of the fit to the V+A 
spectral moments reads 

as(TO2) = 0.344 ±0.005 ±0.007 , (9) 

where the first error is experimental and the sec- 
ond is theoretical. When evolving this value to 
the Z scale P(see Fig. [2|) one gets 

ai^^ (m|) = 0.1212 ± 0.0005 ± 0.0008 ± 0.0005 , 

(10) 

where the first two errors are propagated from ^ , 
and the last one summarises uncertainties in the 
evolution. The consistency between this result 
and the value found by a global fit to elec- 
troweak data at the Z-mass scale |T], a'"j\m'^) — 

ai^\ml) = 0.0021 ± 0.0029, provides the most 
powerful present test of the evolution of the 
strong interaction coupling, over a range of s 
spanning more than three orders of magnitude. 

5. CONCLUSIONS 

Motivated by some new results both on theoret- 
ical and experimental grounds, we have revisited 
the determination of asim^) from the ALEPH 
T spectral functions. We have re-examined two 
common numerical methods: we have identified 
specific consistency problems of FOPT, which do 
not exist in CIPT. The r measurement of as 
evolved to the Z scale is found to be in excellent 
agreement with the direct determination from Z 
decays. Both results are the only ones at N^LO 
order so far, confirming the running of as between 
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Figure 2. Top: The evolution of Q;s(m^) to higher 
scales /i using the four-loop RGE and the three- 
loop matching conditions applied at the heavy 
quark-pair thresholds (hence the discontinuities 
at 2rnc and 2rnii). The evolution is compared 
with independent measurements covering fi scales 
that vary over more than two orders magnitude. 
Bottom: The corresponding as values evolved to 
mz ■ The shaded band displays the r decay result 
within errors. 

1.8 and 91 GeV, as predicted by QCD, with an 
unprecedented precision of 2.4%. 
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